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We compare the solutions to the short baseline neutrino anomaly based on oscillations to sterile 
neutrinos and the Stochastic Neutrino Mixing Mechanism (SNMM) through an analysis of the 
present neutrino data. The SNMM suggests worse fits than a 3 + 1 sterile neutrino model, although 
it cannot be discarded by present data. We propose an experiment to distinguish between both 
solutions, based on placing a ®Li source inside a 5kton-yr detector (like SNO). We studied the 
sensitivity of such an experiment, which makes it possible to discriminate within 2a the SNMM 
from the 3+1 sterile hypothesis for some particular values of the relevant parameters in 5 kton-years 
of running. 


I. INTRODUCTION 

The recent analysis provided by Ref. [T] increased the 
theoretical prediction of the reactor anti-neutrino mean 
flux from ^^®Pu, ^^^Pu, and by about 3 per¬ 

cent. The published reactor experiments at short base¬ 
line (< 100 m) reported a ratio of observed event rate over 
predicted rate of 0.976 ± 0.024. In the new approach, this 
ratio shifts to 0.943 ± 0.023, leading to a deviation from 
unity at 2.5cr. This is called the reactor anti neutrino 
anomaly [5]. The ratio of observed over predicted events 
was also studied with an intense artificial ^^Cr and ^^Ar 
radioactive sources, which were placed inside the detec¬ 
tors of GALLEX [21 m and SAGE [5] , and also indicated 
rates smaller than unity of 0.86 ± 0.05, which has a devi¬ 
ation from the unity at 2.8 a. This is called the Gallium 
anomaly |6]. 

Moreover, searches for ^ Df, oscillations conducted 
by the Liquid Scintillator Neutrino Detector (LSND) 
found an excess of events above the expected back¬ 
ground |7]. The MiniBooNE experiment also reported 
an excess of events at 3.8 cr, combining the Ue and i>e 
data sets |S]. 

This panorama cannot be explained by the standard 
model of neutrino oscillations, in which only active neu¬ 
trinos are allowed to oscillate. Once the oscillation 
length related to the mass eigenstates 1-2 is well mea¬ 
sured by Solar and KamLand |5] data and the oscilla¬ 
tion length of mass eigenstates 1-3 is also well established 
by atmospheric neutrinos |10j and long baseline experi¬ 
ments mm, the excess in MiniBooNE, LSND, old re¬ 
actors and Gallium anomalies, cannot be accommodated 
in a three oscillating neutrino families context. 

To explain such anomalies, the neutrino sterile hypoth¬ 
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esis is the most investigated scenario m- Sterile neutri¬ 
nos are additional states beyond the standard electron, 
muon, and tau flavors, which do not interact by charged 
currents (GG) neither by neutral currents (NG). Sterile 
neutrinos are connected to additional mass states via an 
extended mixing matrix with extra mixing angles. 

Other proposals to handle such anomalies were also 
studied by several authors, like GPT violation m and 
extra dimensions |15] . In special, a proposal was made 
in [16] , based on what we called Stochastic Neutrino Mix¬ 
ing Mechanism (SNMM) in which the mixing angles are 
stochastic variables that can be different in the process 
of neutrino creation and detection. We argued that such 
mechanism seems to supply a possible explanation to all 
these anomalies. 

In this paper we will improve the analysis made in |16j 
adding the MiniBooNE data and making a spectrum 
analysis of the accelerator data. We also will compare 
the SNMM with the sterile hypothesis through a statis¬ 
tical analysis. 

We start presenting the 3+1 sterile neutrinos model 
analysis m- Then we introduce the foundations of the 
SNMM and the dependence of their free parameters with 
the energy. We compare, through a statistical analysis, 
the fits of sterile neutrinos and the SNMM. In the last 
section we study a ®Li source as a possible experimen¬ 
tal test to distinguish the solutions to the short baseline 
anomalies based on the SNMM and sterile neutrino hy¬ 
pothesis. 

II. FRAMEWORK OF STERILE NEUTRINO 
AND STOCHASTIC NEUTRINO MIXING 
MECHANISM 

A. Sterile Neutrinos 

While the indications of sterile neutrino oscillations 
have historically been associated only with appearance- 
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based short baseline experiments [7], the recently obser¬ 
vations of disappearance of r'e and Pg laii!, and new 
data reporting appearance of Pg [H] ) provides further mo¬ 
tivation for these models. The more economical model 
is the so-called 3-1-1 model which involves a fourth 
neutrino states, sterile, which relates to the usual active 
eigenstates by a 4 x 4 mixing matrix, Va = ^ai^i where 
i=l,2,3,4. 

Assuming very-short baselines, where the only relevant 
mass-scale is the one involving the fourth neutrino family, 
the flavor oscillation probabilities in sterile neutrino case 
are given by: 



for the electronic neutrino survival probability, and: 


point to explain the data related to short baseline anoma¬ 
lies. This is similar to non-unitary oscillation mechanism 
described in Ref. m- 

The stochastic nature of our mechanism appears as we 
integrate over the stochastic variables 6c and 9d, which 
are the neutrino mixing angle at the creation and de¬ 
tection respectively. The total oscillation probability be¬ 
comes: 

^7r/2 ^7r/2 

< P{vc -t vp) >= / Pu^^vi}f{9c)f{0d)d9cd9d, 

JO JO 
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where the <> symbols represent the averaging due to 
the stochastic mechanism. Here a and /3 denotes the 
neutrino flavor e, ^ or r, f{9) is a distribution function 
to be dehned below, and: 

7 
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for the transition probability. Here L is the neutrino 
traveled distance, E is the neutrino energy. Therefore, 
all oscillation data of short baseline experiments can be 
described by adding these parameters: and C /^4 ele¬ 
ments of the mixing matrix and the mass scale . 

Nevertheless, introducing sterile neutrinos can have 
implications in cosmological observations, especially in 
measurements of the radiation density in the early uni¬ 
verse if the extra neutrinos have signihcant mass (>1 
eV) and do not decay. Recently, Ref. [T5] estimates the 
effective number of neutrinos to be = 3.30 ± 0.27, 
that indicates disfavored limits to sterile neutrinos. A 
complete analysis of global fits to sterile neutrinos was 
obtained in m- Despite the fact that these fits pro¬ 
vide excellent results to the anomalies, it is evident that 
the great values associated with these sterile states mass 
are in conflict with cosmological bounds [5D]. Having 
in mind this incompatibility of the sterile neutrino sce¬ 
nario with cosmological data, besides the fact that there 
is no direct evidence of sterile neutrinos, we start study¬ 
ing another proposal which does not add any new massive 
state, therefore not violating cosmological constraints: 
the SNMM. 


B. Stochastic Neutrino Mixing Mechanism 


where 7 and S run from 1 to 3 and is the usual 

squared mass difference between the mass eigenstates in¬ 
volved in the oscillation process. 

In order to parametrize this different mixing angles in 
the creation and detection procedure we define a 3 x 3 
mixing matrix at the moment of the neutrino creation 
{U^) and at the detection moment in the following 

way: 

( c.d c.d 

C 42 C 13 S 42 C 13 S 13 

— Sj 2 '^C 23 — 042*^523813 — Si 2 ^S 23 Sl 3 S23C13 

■*12^®23 — d{’2 C23S13 —dl2 S23 — C23S13 C23C13 


where Cij = cos 9ij , Sij = sin 9ij , the stochastic variables 

c,d /ic,d 1 c,d • /qc.d i /ic,d 

are c {2 = cosu -^2 ^12 ~ ^12 

sume values in the interval [ 0 , 7 r/ 2 ]. For the general case 
in which 9 12 ^ 9^2, the survival probability, Eq. ([^ is 
less than unity even at zero distance, also the conversion 
probability is non-zero for zero distance. Such behavior, 
which is not allowed in the usual oscillation processes, is 
the essence of the SNMM. For simplicity, we choose the 
stochastic variables distribution function as a gaussian 
distribution: 


fi^i 2 ) = exp 



( 5 ) 


The Stochastic Neutrino Mixing Mechanism make two 
assumptions: 

1 . the mixing angles that compose the neutrino eigen¬ 
states at creation and detection can be different; 

2 . the mixing angles are stochastic variables. 

The first assumption can generate a non-zero oscillation 
probability even for zero distance, which is the central 


where the N is the normalization factor. We assume for 
a an energy dependence given by a = A -|- B with 
A, B and n as free parameters. 

In order to clarify the model we can write explicitly 
the conversion probabilities replacing the coefficients: 

^a,7(5 —^ Uctd ^ 

= J d9'i2f{0l2)K^Ks = J d9UiOM^U^s ( 6 ) 
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in Eqs. ^ and ([^, obtaining: 

< P(^U(^ ^ ^ ^ ^ ^a,77^/3,77 

7 

/ \ 

+ 2 ^ ^ Va,'YS'V/3,'YS COS I Z/ | , (7) 

J>s \ ) 

where we see that the stochastic effect works as a mod¬ 
ulation in the probability even in very short-baselines 
distances. The second term in Eq. Q modulates the 
amplitude of oscillation, but not destroying the pattern 
of oscillation which is important to describe the acceler¬ 
ator, reactor and the atmospheric data. 

In order to perceive the changes on flavor conversion in¬ 
duced by the stochastic mechanism, is particularly clear 
to analyze the electronic survival probability at very 
short baselines in the approximation @13 = 0. Making 
L —>■ 0 in Eq. Q, under this approximation, we can 
rewrite the survival probability as: 

< P{y^ i/e) >= -f 2 (8) 

= 1 -I- 2 [< C 12 S 12 >^ — < X ®12 >] 

where the <> symbols represent the averaging due to the 
stochastic mechanism, like in Eq. §. From this expres¬ 
sion is easy to see that if the stochastic mechanism is not 
present, due to unitarity the survival probability at very 
small baseline is 1. So, the effect of stochastic neutrinos 
is a zero distance effect, allowing flavor conversions for 
very short baselines. 

For the solar neutrinos, the averaging over oscillation 
takes place, and only the first term in Eq. Q provides 
flavor conversion. And for this term we should analyze 
the stochastic mechanism effects in two different energy 
regimes. For high energy solar neutrinos, the matter ef¬ 
fect dominates over the oscillation term Ato^/ 4E, and 
all neutrinos are produced as mass eigenstate v^. As a 
result, the averaging over the stochastic angle is not ef¬ 
fective and the only effect of the stochastic mechanism 
would be an averaging over a mixing angle on detection 
point, leading to a negligible effect on probability. For 
low energy solar neutrinos, the matter effect are not large 
enough to change the mixing angles, and we obtain the 
first term of Eq. 0 averaged over production and de¬ 
tection mixing angles. We numerically checked that the 
effect of the stochastic mechanism leads to a decrease in 
2 % on solar neutrino survival probability, which would 
only marginally change the data analysis involving these 
neutrinos. For these reasons, we will focus our analysis 
just in experiments of very short baselines. 

III. THE ANALYSIS 

We present the comparison of our proposal, the SNMM 
mechanism, with the sterile neutrino mechanism and 
the usual three neutrino scenario for the following data: 



Su 

3-bl 

SNMM 

data points 

Reactors 

34.41 

22.58 

30.16 

21 

SAGE /GALLEX 

8.09 

5.26 

3.27 

4 

LSND i? 

16.48 

3.77 

3.89 

8 

MiniBooNE i/ 

18.69 

6.98 

15.54 

11 

MiniBooNE 1 / 

28.56 

11.76 

20.81 

11 

NuTeV P 

25.32 

25.32 

25.84 

17 

NuTeV ly 

16.59 

16.58 

10.18 

17 

Total 

148.14 

92.24 

109.72 

89 

d.o.f. 

89 

86 

86 

- 

G.O.F. 

8.4 X IQ-'^ 

0.303 

0.043 

- 


TABLE I: best fit values for each experiment for the 

combined sets for the 3-1-1 sterile model, for the stochas¬ 
tic mechanism and for the usual iv case. We also show 
the degrees of freedom (d.o.f) and the goodness of fit 
(G.O.F.) for the combined data. 

all set of old reactors (21 data point) [2], SAGE and 
GALLEX calibration experiments (4 data point) [3H5], 
anti neutrino channel in LSND (8 data point) [7], neu¬ 
trino and anti neutrino channel in MiniBooNE (22 data 
point) [8] and neutrino and anti neutrino channels in 
NuTeV (34 data point) [55], using a defined as: 

89 

yX^(i?*-i?f)W,-i(i?*-i?|), (9) 

i=l 

where the i?* are the theoretical predictions, i?® are 
the experimental measurements, W is the correlation 
function, and the sum is performed using 89 data 
points in total. For the old-reactors data and for the 
GALLEX/SAGE data we use off-diagonal elements of the 
correlation matrix from Ref. [Tj. 


A. Comparing Sterile and SNMM 

We show our y^ analysis results for all data and three 
different scenarios: the 3-1-1 sterile model, the Stochastic 
Neutrino Mixing Mechanism and the three-neutrino case 
in Table |T] We present the best fit values for y^ for each 
experiment and also for the full data set. 

For the three neutrino case, one cannot have short 
baseline oscillation using the standard values of mixing 
angles and mass differences for any of the experiments 
that we analyzed. For the experiments that show no 
oscillation, such as NuTeV, there is a good agreement 
between data and theoretical predictions. But for other 
experiments, such as the reactors and the Mini-BooNE 
experiment, there is a disagreement. The y^ has no free 
parameters and gives a very bad description of the data 
with the value of 148.14 with a very bad goodness-of-fit 
of 8.4 X 10"^ 

The sterile neutrino parameters that produce the best 
fit to data for the 3-1-1 model are Am^i = 0.42 eV^, 
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FIG. 1: The black line represents the behavior of the 
3+1 sterile neutrino model using the best fit values ob¬ 
tained in the global fit (Amli = 0.42 eV^, U;j 4 = 0.29 
and Ue 4 = 0.14), the red line represents the behav¬ 
ior of the standard Sv scheme using the best fit values 
of |23| and the blue line represents the behavior of the 
SNMM using the best fit values obtained in the global 
fit (a = A + B/E", Ab.f. = 0.026, Bb.f. = 0.26 and 
rib.f. = 0.80). All of the results are obtained using a fixed 
distance of 10 m. P 3+1 depends on LjE, and is straight¬ 
forward to re-escalate the figure for different values of L, 
while the Psnmm has the same dependence with L as 
P 3 ., which is negligible for very-short baselines. 


U ^4 = 0.29 and Ue 4 = 0.14. Our values are in reason¬ 
able agreement with the values of other 3+1 model analy¬ 
ses [13]. The goodness-of-ht value shows that the sterile 
neutrino case is a reasonable explanation for the com¬ 
bined data with the best fit value providing Xhf =92.24 
for 86 d.o.f, with a goodness-of-ht value of 0.303. 

For the Stochastic Neutrino Mixing Mechanism mech¬ 
anism, the free parameters are the energy dependence 
of the gaussian width a = A + B/E'^. We have found 
for the global analysis that Ab.f. = 0.026, Pb.f. = 0.26 
and Uh.f. = 0.80 for neutrino energies given in MeV. The 
SNMM ht is better than the 3z/ case, but it is not bet¬ 
ter than the sterile case. The main difference in be¬ 
tween these two scenarios comes from the reactor anal¬ 
ysis, mainly because the flat pattern reproduced by the 
SNMM does not ht so well the data as the rippled pat¬ 
tern of the sterile neutrino model. The same behavior 
also provides better results to the sterile neutrino in the 
region of MiniBooNE data. In order to clarify this we 
show in Fig Q the behavior of the 3iy, the sterile and 
SNMM cases for the as a function of the energy 

for length of 10m, which is a characteristic length for re¬ 
actors experiments. In this hgure the standard 3 v model 
predicts none oscillation while the SNMM provide a con¬ 
version effect and the sterile neutrino have an oscillation 
effect. For the best fit values of sterile case and SNMM 


we have a different conversion pattern, the sterile case 
have a short oscillation length and the SNMM have a 
conversion effect due the averaging effect in Eq. 0. 


IV. POSSIBLE TESTS 

One of the possible tests to confirm or exclude the 
SNMM in comparison to the sterile neutrino models con¬ 
sists in allocating a source of ®Li inside a detector sensi¬ 
tive to charged and neutral current signals, such as SNO 
experiment [24] . for example. The ®Li source produces 
i>e neutrinos with the average energy of 6.4 MeV and we 
will assume a detector located at most at 16m from the 
source. Using the SNO setup that consists in a tank filled 
with heavy water, reactions like + D ^ n + n + e^ and 
D + D —>■ i/' + n+p, will happen. For these configuration, 
with LjE ^ 16/6.4 ~ 2.5m/MeV, the standard oscilla¬ 
tion scenario predicts no-oscillation and any oscillation 
seen is due to new physics. 

The original motivation for this experiment is to ob¬ 
serve the oscillation pattern of the sterile neutrino model 
for the typical parameters found in previous section. 
We suggest to extend the range of possible models to 
be tested in this experiment, including the Stochastic 
Neutrino Mixing Mechanism, by a measurement of the 
charged/neutral currents of Vf.. 

For this we computed the expected rate for the stan¬ 
dard 2>v scenario both for the charged current, as 

well for the neutral current as 

„cc/NO ^ I (10) 

where the cross section (ccc/nc) is obtained in [25j and 
the spectrum (S'sLi) is provided by [26]. Here nr is the 
number of targets and T is the lifetime. 

For the sterile mass model, we expect to see oscilla¬ 
tions, and then the number of events for charged and 
neutral currents is modified as 




UtTSsu ( (Tec 
IP- I CTnc 



dEdL. 


Here 'Ylix ^ 1 presence of sterile neutrino. 

For the Stochastic Neutrino Mixing Mechanism we 
have similar expressions. 


( ^SNMM 

V ^SNMM 


utTSsi^i / (Tec 
L2 \ (TnC 


pSNMM \ 

j dEdL, 


where in SNMM we have = 1. In the neu- 

trino sterile hypothesis a detection rate decrease is ex¬ 
pected in both neutral current and charged current chan¬ 
nels due to the oscillation of electronic anti-neutrino in 
sterile anti-neutrinos. But in the SNMM hypothesis, 
while the charged current detection rate decrease due to 
the SNMM zero distance effect, the total number of ac¬ 
tive neutrinos remains constant, and then the detection 
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FIG. 2: Ratios (f)i and (pi defined in Eq. (111. To plot 


the 3+1 curves was used the best fit values of our sterile 
analysis. For the Stochastic Neutrino Mixing Mecha¬ 
nism we use the la values of A,B and n parameters of 
our parametrization of Gaussian width a = A-\- BjE^. 
The points indicated by amax,ab.f. and amin indicate 
the maximum, best ht values and minimum values for 
Stochastic Neutrino Mixing Mechanism. 


rate through neutral current remains the same as in the 
standard Siz oscillation model. 

To characterize these differences between the sterile 
model, the Stochastic Neutrino Mixing Mechanism and 
the standard 3 jz model we propose three observables: The 
ratio of observed charged current events, and the 

neutral current events, to the standard 3 jz case and 
the double ratio of NC over CC ratios are defined as, 
respectively, 


4>i 


Ng^' 




<^2 
</>l ' 


( 11 ) 


We will assume an experiment running for 5 years with 
the detector of the size of SNO detector, 1 kton. We 
can do two types of analysis, the rate and the shape of 
measured data for this configuration. 

We show in Fig. |^the predicted values for cpi and (j )4 
for the best fit values of 3+1 model and for the Stochastic 
Neutrino Mixing Mechanism. We can notice that from 
the rate only analysis we cannot discriminate between 
the standard 3;/ scenario and the best fit of 3+1 sterile 
model since the predictions overlap within Icr. For the 
Stochastic Neutrino Mixing Mechanism, we have the Icr 
values of free parameters a = A + BJE'^ of this mech¬ 
anism to be Ajnin = 0.02, Amax = 0.028, Bmin = 0.26, 
Bmax = 0.36 and rimin = 0.78 and Umax = 0.79. One can 
see that we can discriminate at 2a the standard case 
from the Stochastic Neutrino Mixing Mechanism. Even 
more important, one can see from Fig. ([^ that for some 
values of the a parameter close to Omax, it is possible to 


distinguish SNMM from the 3+1 sterile model at 2cr in 
5kton-yr of running. 

This shape analysis was already performed by |26j 
in a KamLand-like detector which has almost the same 
fiducial mass than SNO. Ref. [55] concluded that in a 
five-year running experiment it is possible to distinguish 
between the 3+1 and the standard oscillation model in 
a KamLand detector size. Since the Stochastic Neutrino 
Mixing Mechanism will provide no spectrum shape 
changes, this test will not add any information. A com¬ 
bined analysis (shape + rate) using CC and NC in the 
SNO experiment can definitively distinguish between the 
sterile neutrino model, the standard oscillation model or 
the SNMM. As a disclaimer we use only statistical errors 
and the existence of systematic errors can change the 
conclusions. The rate (pA is less sensitive to systematic 
errors and give a better perspective to discriminate 
between the sterile neutrino and the stochastic models. 


V. CONCLUSIONS 

The Stochastic Neutrino Mixing Mechanism supplies 
an alternative explanation for the short baseline experi¬ 
ments that can explain together the positive oscillation 
signal of reactors, SAGE/GALLEX experiments, Mini- 
BooNE and LSND experiments and the negative results 
of NuTeV. The agreement with data improve compared 
with standard 3iz neutrino scenario, but not with the 
same quality of 3+1 model. 

To distinguish the solution offered by the SNMM from 
that one originated from the sterile neutrino hypothesis 
we can use an artificial radioactive source and monitor 
the rate and shape of the spectrum of produced neutri¬ 
nos. We propose the following variables, the GG ratio, 
NCC jyNC 

the NG ratio (j )2 = and the NC/GG 

02 

double ratio 04 = — to discriminate between the two 
01 

scenarios. 

In an experiment with production of electronic anti¬ 
neutrinos, if a decrease in the NG signal is observed, the 
SNMM will be excluded, otherwise, an evidence in favor 
of the SNMM will be found. A combined rate and shape 
analysis using an SNO-like detector can point out in favor 
of the SNMM or the Sterile Neutrino hypothesis in a 
5kton-yr running experiment. 
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